REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  i 


Form  Approved 
OMB  No.  0704 


*ta  needed  a^icl  complefing'and  revfo^^  Wstarttones!^^  *0Urces>  9atherin9  and  maintaining  the 

this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  ReDorts  (0704^lR^5irt  '2?Tat'0n»lnc,udln9  su99estions  fo,  reducing 

4^?«Ftfpondents  should  awar®  that  notwithstanding  any  other  provision  of  law  no  person  shall  be  subject  to  anv  ne™Hu  L-  1 21  ?  JeJfrson  Davis  Highway,  Suite  1 204,  Arlington,  VA  22202- 

valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  MOVE  ADDRESS  ”  to  any  penalty  for  failing  to  comply  with  a  collecbon  of  information  if  it  does  not  display  a  currently 

1.  REPORT  DATE  (DD-MM-YYYY)  I  - -  - - - - - - - - 


30-3-2005 


RF.PRTNT 


4.  TITLE  AND  SUBTITLE 

Monopole-Quadrupole  Model  of  Spacecraft  Charging  in  Sunlight 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S)  “ 

M.  Tautz*  and  S.  T.  Lai 


5d.  PROJECT  NUMBER 

5021 


5e.  TASK  NUMBER 

RS 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratorv/VSBX 

29  Randolph  Road 

Hanscom  AFB  MA  01731-3010 


5f.  WORK  UNIT  NUMBER 

A1 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


AFRL-VS-HA-iTR-2  005- 1003 


.  SPONSORING  1  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/VSBX 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  ~ - 

Approved  for  Public  Release;  Distribution  Unlimited. 

*AER,  Inc.,  Lexington,  MA 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


13.  SUPPLEMENTARY  NOTES 

REPRINTED  FROM  PROCEEDINGS: 


43rd  AIAA  AEROSPACE  SCIENCES  MEETING  AND  EXHIBljr 
10-12  Jan  2005,  RENO,  NV  j 


Geophysics  research  spacecraft  are  often  covered  partially  with  thin  layers  of 
dielectrics  such  as  kapton  or  thermal  blanket  materials.  Analytical  modeling  of  the 
charging  of  a  dielectric  spacecraft  in  sunlight  needs  to  go  beyond  the  ihonopole 
potential  and  include  a  dipole  term.  If  the  capacitances  of  the  dielectric  layers  are 
sufficiently  large  and  the  spacecraft  spin  rate  is  sufficiently  fast  the  potentials  may 
be  determined  on  the  basis  of  spin  averaged  spectra  of  electron  and  ion  fjfuxes.  In 
such  a  condition  the  potential  distribution  resembles  that  of  a  monopole- 
quadrupole  system.  We  formulate  the  monopole-quadrupole  case  and  compare  it 
to  the  monopole-dipole  model.  1 


15.  SUBJECT  TERMS  - - - - 

Spacecraft  charing  Monopole-dipole  model 

Spinning  spacecraft  Charging  in  sunlight 

1$.  SECURITY  CLASSIFICATION  OF: 


Monopole-quadrupole  model 


a.  REPORT 

UNCLAS 


UNCLAS 


17.  LIMITATION 

18.  NUMBER 

OF  ABSTRACT 

OF  PAGES 

c.  THIS  PAGE 

SAR 

UNCLAS 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Shu-T.  Lai 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

781-377-2932 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


AIAA-2005-0068 


AFRL-VS-HA-TR-2005-1 003 

Monopole-Quadrupole  Model  of  Spacecraft  Charging  in 

Sunlight 


Maurice  Tautz# 

AER,  Inc.,  Lexington,  MA  01 730 

and 

Shu  T.  Lai  ** 

Space  Vehicles  Directorate,  Air  Force  Research  Laboratory,  Hanscom  AFB,  MA  01 731 


Geophysics  research  spacecraft  are  often  covered  partially  with  thin  layers  of 
dielectrics  such  as  kapton  or  thermal  blanket  materials.  Analytical  modeling  of  the 
charging  of  a  dielectric  spacecraft  in  sunlight  needs  to  go  beyond  the  monbpole 
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Nomenclature 

quadrupole/dipole  strength  relative  to  the  monopole 
coefficients  giving  the  strength  relative  to  the  monopole 


=  monopole  potential 
=  nth  order  Legendre  polynomial 
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I.  Introduction 

Daylight  charging  of  a  spherical  object  covered  with  dielectric  was  first  studied  by  Mandell 
et.al.1,  Higgins2  and  Besse  and  Rubin3.  The  Laplacian  potential  distribution  around  the  satellite  is 
represented  by  a  combination  of  monopole  and  dipole  terms.  Axial  symmetry  around  the  sun  direction 
is  assumed,  so  the  object  could  be  either  stationary  or  rotating.  The  monopole  and  dipole  terms 
combine  to  form  a  potential  barrier  outside  the  sunlit  surface  which  acts  to  suppress  the  escape  of 
photoelectrons,  leading  ultimately  to  current  balance.  In  this  paper  we  describe  a  similar  model,  where 
the  dipole  term  is  replaced  by  a  quadrupole  contribution.  In  the  monopole-quadrupole  system,  the 
satellite  spins  about  the  polar  axes  of  the  sphere  and  the  sun  is  interpreted  to  be  shining  on  the  belly- 
band.  Axial  symmetry  around  the  spin  pole  is  assumed  and  thus  the  satellite  must  be  rotating  rapidly 
so  that  only  revolution  time  averaged  photoemission  is  experienced.  An  axially  symmetric  potential 
barrrier  forms  at  the  belly-band  to  block  photoelectrons,  enabling  the  net  current  to  go  to  zero.  The 
potentials  give  the  first  order  charging  reponse  in  daylight  of  a  roughly  spherical,  non-conducting,  fast 
spinning  spacecraft  in  a  low  density  space  environment. 

To  facilitate  comparison,  the  monopole-quadrupole  and  monopole-dipole  analytic  models  are 
developed  in  parallel.  There  are  two  free  parameters  in  the  models.  The  first  parameter,  K,  is  the 
monopole  potential.  The  second  parameter,  A,  is  the  quadrupole/dipole  strength  relative  to  the 
monopole.  The  basic  equations  of  the  models  and  the  parameters  are  described  in  section  2.  In  section 
3,  numerical  solutions  are  outlined  and  results  are  given.  Section  4  contains  summary  remarks. 

n.  The  Monopole/dipole/quadrupole  Expansion 

The  meaning  of  rapid  spin  is  relative  to  satellite  surface  charging  times.  If  the  spin  period  is  long 
compared  to  the  differential  charging  time  of  dielectric  surface  elements,  the  surface  is  able  to  respond 
to  the  sun  and  the  motion  may  be  considered  slow.  If  the  spin  period  is  short  compared  to  the  time  it 
takes  a  dielectric  surface  to  charge,  the  motion  is  considered  to  be  fast.  In  this  case  the  surface  would 
respond  only  to  spin  averaged  solar  illumination. 

Consider  a  dielectric-covered  spherical  satellite  that  is  rotating  rapidly  in  sunlight,  so  that  only 
spin  averaged  effects  are  important.  If  the  ambient  space  charge  density  is  low,  which  can  be  the  case 
at  geosynchronous  altitudes,  the  potentials  outside  the  satellite  would  be  given  approximately  by  an 
axially  symmetric  solution  to  Laplace's  equation.  In  spherical  coordinates,  the  Laplacian  potentials  can 
be  written  in  the  form  (see  Schwartz4): 


K^AnPn{t) 


r  „=i 


n= 0 


(1) 


where  r  is  the  radius  and  t  is  the  polar  angle.  The  sum  is  over  n  =  0,  1,  2...  and  Pn(  t  )  is  the  n,h  order 
Legendre  polynomial.  The  constant  coefficient  K  is  the  monopole  potential  and  the  coefficients  An 
give  the  strength  relative  to  the  monopole.  If  we  keep  the  three  lowest  order  terms,  the  potential  will 
be 


V(r,t)  = 


K 


^  |  AW)  !  A2P2{t') 


(2) 


where  A,  is  the  dipole  strength  and  A2  represents  the  quadrupole.  If  we  assume  a  unit  sphere,  we  have 
on  the  surface 


VM=K(\+Am+Am)  <3) 

The  monopole-dipole  case  corresponds  to  Ai=  -A,  A2= 0  and  the  monopole-quadrupole  case  is 
specified  by  A\=  0,  A2  =  A.  We  note  that  the  minus  sign  in  the  monopole-dipole  case  is  arbritary  and 
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agrees  with  previous  treatments:  a  plus  sign  would  put  the  sun  at  the  opposite  spin  pole.  The  K  and  A 
parameters  depend  on  the  balance  of  the  incoming  and  outgoing  satellite  surface  currents  and,  since 
the  ambient  currents  are  here  not  assumed  to  be  known,  are  free  parameters  of  the  models. 

The  overall  shape  of  the  potential  distributions  can  be  seen  by  looking  at  the  north,  middle  and 
south  surface  potentials  at  the  polar  angles  t  =  0,  90  and  180  degrees  respectively.  These  potentials  are 
shown  below,  with  the  monopole-quadrupole  system  on  the  left  and  the  monopole-dipole  case  on  the 
right,  enclosed  in  curly  brackets 


F(1,0)  =  VN  =  K{\+  A) 

{  K(l-A)} 

(4) 

V(\,90 )  =  VM  =  K  (1  -  A  12) 

{  K  } 

(5) 

F(l,180)  =  Vs  =  K{  1  +A) 

{  K(\+A) } 

(6) 

In  these  equations,  we  normally  have  K  <  0  for  negative  voltage  charging  and  A  >  0  for  sunlit 
charging.  The  monopole-quadrupole  system  has  equal  potentials  at  the  poles,  and  a  lower  (negative) 
potential  at  the  belly-band.  In  contrast,  for  the  monopole-dipole  system,  the  potentials  increase 
(negatively)  in  going  from  the  north  to  south  poles. 

The  lower  limit  on  A  is  determined  by  the  charging  threshold,  discussed  in  the  next  section.  A 
practical  upper  limit  on  A  can  be  taken  as  the  point  where  the  lowest  potentials  go  to  zero.  This  gives 
A  =  1  for  the  monopole-dipole  system  and  A  =  2  for  the  monopole-quadrupole  case. 

III.  Solution  of  the  Models 

Photoemission  currents  (positive)  usually  dominate  ambient  currents  in  the  magnetosphere,  yet 
negative  charging  is  reported  in  daylight.  A  mechanism  for  sunlight  charging  is  well  known.  The 
shaded  surfaces  charge  up  and  the  fields  wrap  around  the  object  to  the  sunlit  side.  A  potential  barrier 
forms  outside  of  the  sunlit  surface  to  trap  escaping  photoelectrons,  allowing  current  balance.  The 
models  are  based  on  this  barrier  dominated  scenario.  To  obtain  an  analytic  solution  to  the  problem  we 
neglect  a  self-consistent  photosheath,  which  would  require  Poisson's  equation  and  particle  tracking. 

For  a  potential  barrier  to  form,  we  have  the  condition 

dV(r,t)/dr  =  0  (7) 

and  using  the  previous  expression  (2)  for  V(  r,  t ),  we  get 

4f,-MW)+3M(£)U  (8) 

r  \  r  r  ) 

Solving  this  equation  for  r  gives  the  barrier  radius  Rh(  A,  t ).  Here  and  below  we  give  the  result  for  the 
monopole-quadrupole  system  (the  corresponding  quantities  for  the  monopole-dipole  model  are 
summarized  in  Table  1) 


R,(< 4,t)  =  (-HA)"2 


A  maximum  barrier  radius,  RB,  occurs  at  some  angle  tB.  For  the  monopole-quadrupole  model,  Rb  is  a 
maximum  when  -P2  is  at  a  maximum,  which  occurs  at  tB  =  90°  and  we  have  Rn  -  (3/2  A  )m  .  We 
can  interpret  the  angle  tB  as  specifying  the  sun  direction.  A  barrier  forms  outside  the  sphere  for  ,4- 
values  above  a  threshold  which  can  be  determined  using  the  condition  RB  =  1,  which  gives  A  =  2/3. 
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Figure  1.  Maximum  barrier  radius  as  a 
function  of  strength  A  for  the  dipole  and 
quadrupole  models.  The  barrier  radius 
increases  monotonically  with  A  from  the 
threshold  value  of  unity. 


Figure  2.  Maximum  barrier  height  VB 
(normalized  by  K)  for  the  dipole  and  quadru¬ 
pole  models.  The  barrier  heights  monotonically 
increase  from  zero  at  the  thresholds. 


Figure  1  shows  the  maximum  barrier  radius  for  the  two  models  in  the  range  of  interest.  The  barrier 
radius  increases  monotonically  with  A,  from  the  threshold  value  of  one.  The  maximum  height  of  the 
potential  barrier  is  given  by 


VB=V(JtB,tB)-V(l,tB) 

For  the  monopole-quadrupole  case,  we  find  a  normalized  barrier  height 


Vo\3/2 


Vb=K 


2 

v3y 


A-V2+~- 1 
2 


(10) 


(11) 


A  plot  of  Vq!  K  for  the  two  models  is  given  in  Figure  2.  The  barriers  monotonically  increase  from  zero 
at  the  thresholds. 

We  can  also  determine  the  angular  width,  tw,  of  the  barrier  region,  given  an  -value  above 
threshold.  The  condition  RB  =  1  gives  the  angle  fi  at  which  the  barrier  returns  to  the  surface  and  hence 
tw  =  ±  ( tB  -  ti  ).  Solving  for  the  monopole-quadrupole  case,  we  get 


t]  =  cos 


A- 2/3 


3  A 


(12) 


At  threshold  tw  =  0  and  at  maximum  A  we  have  tw  =  ±  28.1°.  Since  the  maximum  barrier  width  is 
less  than  the  nominal  angular  width  of  photoemission  exposure  ±90°,  we  do  not  have  self-consistant 
photoemission  dynamics. 

The  sun  to  shade  potential  ratio,  Vss,  in  the  monopole-quadrupole  case  is  given  by 


_VM  _\-AI2 

K  no  —  - 

vs  1 +  A 


(13) 
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STRENGTH  A 


Figure  3.  The  sun  to  shade  potential  ratio 
as  a  function  of  strength  A  for  the  dipole  and 
quadrupole  models. 

which  at  threshold  is  2  /  5.  A  compari-son 
plot  of  Vss  is  given  in  Figure  3.  Since  Vss 
for  both  cases  is  monotonically  decreasing, 
the  models  predict  that  this  ratio  will  be 
below  the  threshold  values. 

Figure  4  shows  the  monopole-quadrupole 
potentials,  normalized  to  K,  in  the  physical 
space  surrounding  the  sphere.  The  contour  plot 
shows  a  Y  =  0  slice  of  data  expressed  in  X,  Y,  Z 
coordinates,  which  are  normalized  to  the  sphere 
radius.  The  A  parameter  has  been  selected  so 
that  the  barrier  potential  ratio  VB  /  K  ~  0.01. 
This  sets  up  a  barrier  of  -10  V,  given  a 
monopole  potential  of  -1000  V.  An  axially 
symmetric  potential  barrier  forms  at  the  belly- 
band.  The  contour  plot  for  the  corresponding 
monopole-dipole  case  is  given  in  Figure  5. 


Figure  4.  Potentials  for  the  Monopole- 
Quadrupole  system  in  the  plane  Y=0. 

The  parameters  used  are  tB  =  90°,  A  =  0.811, 
Rb=1.103,  Vb/K=0.01. 


Figure  5.  Potentials  for  the  Monopole-Dipole 

system  in  the  plane  Y=0.  The  parameter  used 
are  tB=0,  A  =  0.576,  RB  =  1.152,  Vb/K=0.01. 


IV.  Summary 

For  a  spherical  fast-spinning  dielectric-covered  satellite  that  is  located  in  a  low  density  ambient 
space  plasma,  we  have  developed  a  simple  analytic  model,  the  monopole-quadrupole,  for  charging  in 
sunlight.  The  model  sets  up  an  axially  symmetric  potential  barrier  at  the  belly-band  of  the  satellite 
which  acts  to  suppress  photoemission  and  leads  to  current  balance.  We  can  interpret  the  model  as 
representing  spacecraft  charging  with  the  sun  at  right  angles  to  the  spin  axis.  The  belly-band  charges 
less  (negatively)  than  the  spin  poles  and  the  sun  to  shade  potential  ratio  lies  below  its  threshold  value. 

The  monopole-quadrupole  and  the  monopole-dipole  systems  are  similar  in  that  they  ignore 
many  details  of  satellite  construction  and  charging  environment  in  order  to  capture  the  main  effects. 
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For  both,  the  analysis  is  based  on  axial  symmetry  and  Laplacian  potentials.  The  major  differences 
between  the  models  are:  (1)  The  sun  locations  are  90  degrees  apart  (the  belly-band  versus  the  spin 
pole),  (2)  The  charging  thresholds  are  slightly  different  (  A=  2/3  versus  1/2  ),  (3)  The  sun  to  shade 
ratios  are  similar,  but  not  the  same  ( Vss  =  2/5  versus  1/3). 

In  the  monopole-dipole  system  the  dark  side  is  constantly  hidden  from  the  sun  during  a  spin 
revolution  and  this  model  includes  the  slow  spin  limit.  For  the  monopole-quadrupole  system  the 
shaded  side  is  instantaneously  changing  and  the  results  are  only  valid  for  fast  rotation  rates,  where 
spin  averaged  photoemission  becomes  a  good  approximation. 
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